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Abstract 



We compute the leading quark mass dependence of the matrix elements of 
isovector twist-2 operators between octet baryon states in quenched QCD 
using quenched chiral perturbation theory. There are contributions of the 
form m q logm g and m q , in analogy with QCD, but there are also contributions 
of the form logm g from hairpin interactions. The nucleon does not receive 
such hairpin contributions. 
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I. INTRODUCTION 



The parton distribution functions (PDF's) of the nucleon are fundamental quantities 
associated with the strong interactions. Extensive experimental investigations have been 
undertaken during the past three decades to measure these distributions via deep-inelastic 
scattering (DIS) of leptons from protons and light nuclei. Due to the intrinsically non- 
perturbative nature of the strong interactions in the low-momentum region, theoretical 
efforts to understand these distributions have had only limited success. With the ever 
increasing power of computers and significant developments in algorithms used to numer- 
ically simulate QCD on the lattice, it is hoped that properties of the PDF's can be de- 
termined from first principle at some point in the not so distant future (for a review 
see Ref. P|). It is the forward matrix elements of twist-2 operators that are computed 
numerically and these matrix elements are directly related to moments of the PDF's. Of 
course, only matrix elements computed with unquenched QCD with the physical values of 
the quark masses, m q , are to be directly compared with experimental data, but at this point 
in time such computations are not possible. All present computations are performed with 
lattice quark masses, m^ tt ', that give a pion mass of rn 1 ^' ~ 500 MeV, and most simula- 
tions are quenched. Despite the fact that quenched computations require significantly less 
computer time, they cannot, unfortunately, be connected to QCD in any way. While solid 
progress is being made toward unquenched calculations |2| of the moments of the PDF's, 
it is likely that partially-quenched || computations will first provide a reliable connection 
between lattice computations and nature by allowing for calculations with smaller m q and 
thereby minimizing the impact of the m^-extrapolation. However, it is of interest, from 
a theoretical standpoint, to know the moments of the PDF's in quenched QCD (QQCD). 
To extrapolate from m 1 ^' down to m q , the m^-dependence of the matrix elements is re- 
quired ||, and recently chiral perturbation theory (%PT) has been used to determine the 
leading m^-dependence in an expansion about the chiral limit in QCD ||[^] and large-iV c 
QCD |J. %PT also allows for a connection to be made with existing convolution models of 
PDF's, and in addition, shows how to make them consistent with QCD ||. In this work we 
determine the leading m^-dependence of the matrix elements of isovector twist-2 operators 
about the chiral limit in QQCD using quenched chiral perturbation theory (Q%PT) |T0|-[14 



II. Q X PT 

The lagrange density of QQCD is 

£= £ t \i$-m q \ h a q b + £ f [ilp-m q \\qr h 

a,b=u,d,s a~b=u,d,s 

= J2 Q j [ty-no]) Qk , (1) 

j,k=u,d,s,u,d,s 

where q are the three light-quarks, u, d, and s, and q are three light bosonic quarks u, d, and 
s. The super-quark field, Qj, is a six-component column vector with the three light-quarks, 
u, d, and s, in the upper three entries and the three ghost-light-quarks, u, d, and s, in the 
lower three entries. The graded equal-time commutation relations for two fields is 
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Q?(x)Qf(y) 



Qf(y)Qf(x) = rt^ 3 (x-y) 



(2) 



where a, (3 are spin-indices and i,j are flavor indices. The objects rjk correspond to the 
parity of the component of Qk, with rjk = +1 for k = 1, 2, 3 and rjk = for k = 4, 5, 6. 
The diagonal super mass-matrix, rriQ, has entries tuq = diag(m u , m^, m s , m u , m^, m s ), i.e. 
rriu = m u , mg = and m s = m s , so that the contribution to the determinant in the path 
integral from the q's and the g's exactly cancel. 

In the absence of quark masses, the lagrange density in eq. ([!]) has a graded symmetry 
U(3\3)l <g> U(3\3)r, where the left- and right-handed quark fields transform as Ql — > UlQl 
and Qr — > UrQr respectively. However, the functional integral associated with this La- 
grange density does not converge unless the transformations on the left- and right-handed 
fields are related, sdet (Ul) = sdet where sdet() denotes a superdeterminant P|, |TT| , |T2"|] , 
leaving the theory to have a symmetry [577(3 1 3) l g) SU(3\3)r] x U(1)v, where the "x" de- 
notes a semi-direct product as opposed to a direct product, "(g)". It is assumed that this 
symmetry is spontaneously broken [SU(3\3)i <g> SU(3\3)r] x U(1)v — > SU(3\3)v x U(l)v so 
that an identification with QCD can be made. 

The pseudo-Goldstone bosons of QQCD form a 6 x 6 matrix, $, that can be written in 
block form 



X' 
X 7T 



(3) 



where 7r is the 3x3 matrix of pseudo-Goldstone bosons including the rj' with quantum 
numbers of qq pairs, tt is a 3 x 3 matrix of pseudo-Goldstone bosons including the ff with 
quantum numbers of qq pairs, and x is a 3 x 3 matrix of pseudo-Goldstone fermions with 
quantum numbers of qq pairs, 
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K+\ 






Vd 


K° 
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K+ 




Vd 
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^ Xrju 


Xtt+ 


Xk+\ 
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Xtt- 
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{xk~ 




Xrjs ) 



(4) 



As the object 



str ($) 



1 

72 



V 



(5) 



is invariant under [577(3|3)x (g) SU(3\3)r] x U(l)y the most general lagrange density that 
describes low-momentum dynamics will contain arbitrary functions of $o HTl~l , |T2! . At lowest 
order in the chiral expansion, the Lagrange density that describes the dynamics of the 



pseudo-Goldstone bosons is, using the notation of Ref. [TJ 
P 



C 



str 



+ A str [ wlqL + m} Q Y) 1 + a$ 9 M $ ^$ 



Ml §1 



(6) 



where the parameter A is chosen to reproduce the meson masses, and £ is the exponential 
of the $ field, 



£ = exp 



(7) 



3 



With this normalization, / ~ 132 MeV in QCD. In addition, it is understood that the 
operators with coefficients a$ and Mq , the hairpin interactions, are inserted perturbatively. 
Expanding out the Lagrange density in eq. @ to quadratic order in the meson fields, one 
finds relations between the meson masses in the isospin limit, 



1m\ 



m 



i id 



(8) 



The Lagrange density in eq. (§) has been used to compute several observables in the me- 
son sector, such as fx, fn [0 and the meson masses, to one-loop in perturbation the- 
ory []liyr2"[|l6| (for a review see Ref. |T7|1 ). 

The inclusion of the lowest-lying baryons, the octet of spin-1 baryons and the decuplet 
of spin-| baryon resonances, is detailed in Ref. fT3" |. An interpolating field that has non-zero 
overlap with the baryon octet (when the ijk indices are restricted to 1, 2, 3) is fll3 | 



»7 



)Ct,a ^7 ,c ,b 



(9) 



where C is the charge conjugation operator, a,b,c are color indices and a,j3, 7 are spin 



indices. Dropping the spin index, one finds that under the interchange of flavor indices |T3 
-) 1+Wk B ikj , B ijk + H 1+J7<% B jik + (-y+^+w k +mm Bkj . 



B 



:jk 







(10) 



The object B^k describes a 70 dimensional representation of SU (3\3)v [0- It is convenient 
to decompose the irreducible representations of SU(3\3)v into irreducible representations 
SU(3) q (g) SU(3)q ® U(l) 0. The subscript denotes where the "577(3)" acts, either on the 
g's or on the g's. The ground floor of the 70-dimensional representation contains baryons 
that are comprised of three quarks, qqq, and transforms as an (8, 1) of SU(3) q (g> SU(3) q . 
The octet baryons are embedded as 



Babe — -y= ( tabd + eacdBfr ) 



(11) 



where the indices are restricted to take the values a,b,c = 1, 2, 3 only. The octet baryon 



matrix is 



B 



1 vo 
73^ 



V 



J_A - 
^0 



P 
n 

-4rA 



(12) 



The first floor of the 70-dimensional representation contains baryons that are composed of 
two quarks and one ghost-quark, qqq, and transforms as an (6, 3)® (3, 3) of SU(3) q ®SU(3)q. 
The tensor representation a s a b of the (6, 3) multiplet, where a — 1, 2, 3 runs over the q indices 
and a, b = 1, 2, 3 run over the q indices, has baryon assignment 



aSl3 — a.S 31 



a s 12 
1 

V2 



— a s 21 

(6) w +3 



V2 a 



0,^22 



SJ 1 



a$23 



a^32 



V2 



(13) 
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The right superscript denotes the third component of q-isospin, while the left subscript 
denotes the q flavor. The tensor representation gt a of the (3, 3) multiplet, where a — 1,2,3 
runs over the q indices and a = 1,2,3 run over the q indices, has baryon assignment 

A _ (3)^-5 _j3 _ (3)^+1 ,3 _ AO ( U \ 

The aSab and the s t a are uniquely embedded into E^ k (up to field redefinitions), constrained 
by the relations in eq. fllPp: 

B ijk = ]l — i-3S jk for % = 4, 5, 6 and j, k = 1,2,3 
I J_ 

®ijk = -- k-3t a e aij - -j= k - 3 Sij for A; = 4, 5,6 and i,j,cr = 1,2,3 . (15) 

As we are only interested in one-loop contributions to observables with ggg-baryons in the 
asymptotic states, we do not explicitly show the second and third floors of the 70. 
An interpolating field that contains the spin-| decuplet as the ground floor is [13|] 



®ijk = o i-3* CT£ <nfc H = jsSik for j = 4, 5, 6 and i,k,a= 1, 2, 3 



rr-Ot,H 



QTQfQT + Qi' b Q]' c Qk a + QT c Qf a Qt' b ] e abc (CY)„ , (16) 



,j -vie > -°6 k ' -vj 

where the indices i,j, k run from 1 to 6. Neglecting spin indices, one finds that under the 



interchange of flavor indices [13 



% jk = (-) 1+r >^T jik = (-f+^T lkj . (17) 

%j k describes a 38 dimensional representation of S77(3|3)y, which has a ground floor trans- 
forming as (10, 1) under SU(3) q ® SU(3)q with 

%ibc = T a b c , (18) 

where the indices are restricted to take the values a,b,c = 1, 2, 3, and where T a f, c is the 
totally symmetric tensor containing the decuplet of baryon resonances, 

1 1 o 

Till = A , T 112 = j ^122 — ' ^222 = A 

Tn 3 - -^f , T 123 - -/=£ , T 223 - -j=Z 

T 133 = , T 233 = -i=E*'- , T 333 = A" . (19) 

The first floor transforms as a (6, 3) under SU (2>) q ® SU (2>) q and has a tensor representation, 
gXy, with baryon assignment 



The embedding of g,Xij into %j k is unique (up to field redefinitions), constrained by the 
symmetry properties in eq. ([17]): 



7t 

2 ijk 



1 ijk 



+ 



V3 
f 

1 



isXjk for i = 4,5,6 and j, k = 1,2,3 
3Xjfe for j = 4,5,6 and i, A; = 1,2, 3 



^ijjfc = fc-3^ii for fc = 4,5,6 and z, j = 1,2,3 



(21) 



The free Lagrange density for the Bijk and %j k fields is |fT5 , |H |, at leading order in the 
heavy baryon expansion |18| - 



C = i(BvVB) + 2a M (BBM+) + 2(3 M (BM+B) + 2a M (BB) str(7W + ) 

-i (7*vT%) + A (T%) + 2 lM (T»M + %) - 2a M (T%) stx(M+) , (22) 

where M. + = | ^ mg ^t _|_ £mg£), A is the decuplet-octet mass splitting, is the baryon 

four- velocity, and £ = vE- The brackets, (...) denote contraction of lorentz and flavor 
indices as defined in Refs. [jl3| , pT4 . 

The Lagrange density describing the interactions of the baryons with the pseudo- 
Goldstone bosons is [0 

C = 2a (BS^BA^) + 2(3 (BS^A^B) + 2 7 (BS^B) str 

+ 2U{T v S^A i J v ) + f-C[(T v A v B) + (BA,T")] + 2j'(t u S>*%) str(^) , (23) 



where is the covariant spin- vector [p!8|-p0|. Restricting oneself to the qqq sector, it is 
straightforward to show that in QCD 



a = -D + 2F , f3 = ~^D + F 



(24) 



where D and F are constants that multiply the SU(3) q invariants that are commonly used 
in QCD, but it should be stressed that the F and D in QQCD will not have the numerical 
values of those of QCD. In QQCD there is an additional coupling that must be considered, 
7, a hairpin interaction ||13| , that is usually not considered in the xPT description of low- 
energy QCD. In our calculation of the matrix elements of twist-2 operators we will replace 
a and (3 with F and D, but we will keep 7 explicit. In the above discussion, vector and 
axial- vector meson fields have been introduced in direct analogy with QCD. The covariant 
derivative acting on either the B or T fields has the form 

{V»B) ijk = d*B ijk + {V»)\ B ljk + (-)*(»«+*») <y»y? B imk + (-)(vi+vi)(v k +n n ) B . jn > (25 ) 
where the vector and axial- vector meson fields are 



V» = - ( + ) , A" = ~ ( - i^d^i 



(26) 
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III. ISOVECTOR TWIST-2 OPERATORS 



In QCD, the nonsinglet twist-2 operators have the form, 

0& a ...Mn = ^ Xa 70* (< £/*) ... (i S M ) g - traces , (27) 

where the {...} denotes symmetrization on all Lorentz indices, and where A a are Gell-Mann 
matrices acting in flavor-space. The ^ n transform as (8, 1) © (1, 8) under SU(3) L © 

SU(3)r chiral transformations fUTj. Of particular interest to us are the isovector operators 
where A 3 = diag(l, — 1, 0). In the isospin limit, there are no contributions to the matrix 
elements of A 3 from disconnected diagrams. In QQCD the nonsinglet twist-2 operators have 
the form 

Q 0$* ... Mn = i Q T 7{M1 (i B M2 ) ... (i 5 M ) Q - traces , (28) 

where the A a are super Gell-Mann matrices, and we are interested in the isovector ma- 
trix elements with A 3 = diag(l, — 1, 0, 1, — 1, 0) [} The operators Q 0^ 3 2 ^ transform as 
(8, 1, 1, 1) © (1, 8, 1, 1) © (1, 1, 8, 1) © (1, 1, 1, 8) under SU{3) Lq © SU{3)^ © SU(3) Rq © 
SU(3)ftq graded chiral transformations. 



A. Matrix Elements in the Pion 

At leading order in the chiral expansion, matrix elements of the isovector operator 
-O^]] 3 ^ n between meson states are reproduced by operators of the form [|| 

Q 0^... Un - a™ (i) n J - ( — ] str [ StA 3 S m 9 e^ n E + EA 3 5 M1 5 M2 ... 9 Mn Et 

(29) 

where A x = Anf is the scale of chiral symmetry breaking. Operators involving more deriva- 
tives on the meson fields or more insertions of m q make contributions to the matrix elements 
that are higher order in the chiral expansion ||. In QCD, the one-loop diagrams shown in 
fig. |l] with only the octet mesons in the loop give rise to the leading non-analytic terms 
in m q , of the form m g logm g . Such diagrams have the potential to give analogous contri- 
butions to the matrix elements in QQCD. However, explicit computation proves that such 
non- analytic contributions to the matrix element of ^O^'^ ^ n resulting from the diagrams 




1 In principle, the most general diagonal, isovector, supertraceless charge matrix that could be 
considered is A = diag(l, —1, 0, x, y, —(x + y)) (for a discussion about the extension of electro-weak 
operators contributing to nonleptonic decays in QQCD see Refs. [p^j2^| ). For arbitrary x and y 
there will be contributions from disconnected ghost quark diagrams. However, for y = —x there 
will be no contributions from disconnected diagrams involving quarks or ghost-quarks in the isospin 
limit. In this work we have chosen to examine the case where x = — y = +1. 
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Jl - li&'i - JL Jl - 1^1 - JL 
(a) (b) 



It I J % 



(c) 

FIG. 1. One-loop diagrams that contribute to the matrix elements of^O^^, ...^ n in the pion. 
The crossed circle denotes an insertion of an operator from eq. $2f\), arising directly from the twist- 2 
operator. The smaller solid circle denotes an insertion of a leading order strong-interaction vertex 
from eq. ^). Diagrams (a) and (b) are vertex corrections while diagram (c) denotes wavefunction 
renormalization. 



shown in fig. vanish, due to a cancellation between the 'V -loops and "x" -loops. There- 
fore, at one-loop order there are no nonanalytic contributions and the matrix element of 
Q O l )?)t ,, in the pion is 

( Q { ££ -J- = 1 4 a(n) (^) ^ 9m-9*. - traces + - > ( 3 °) 

for n-odd, between an initial pion with isospin index a and momentum q^, and a final pion 
with isospin index j3. The matrix elements vanish for n-even and the ellipses denote terms 
that are analytic in m q , or are higher order in the chiral expansion. 



B. Matrix Elements in the Octet Baryons 

At leading order in the chiral expansion, the matrix elements of ®Oj£s^ ^ are described 

by 

+ 7 (n) v^v^.-.v^ ( r T a ) + v {fllVfl2 ...v^ n _ 2 ( T Mn _ 1 A~' + % n} ) 

— traces , (31) 

where the {...} brackets denote complete symmetrization of the enclosed indices, is the 
four- velocity of the heavy baryon, and 

xi^l^re+en) , (32) 

in complete analogy with the work of Refs. ||[7|]. The index contractions denoted by (...) are 
the same as those in eq. (|23|) defined in Refs. ||T3|JT^] . In general, the coefficients a^ n \ f3^ n \ 7 ( - n * ) 



S 



and a"( n ) are not constrained by symmetries and must be determined from elsewhere. How- 
ever for n — 1 they are fixed by the isospin charge of the baryons to be 



(i) 



+2 







(i) 



-1 , 7 (1) 



a 



(i) 







(33) 



An expression analogous to eq. ([31] ) exists in QCD, and can be obtained from eq. (|3T|) by 
restricting the flavor-indices in each sum to take values i = 1, 2, 3 only. However, one must 
bear in mind that, in general, the values for a^ n \ (3^ n \ r y^ 1 ' and in QQCD will differ from 
those in QCD for n 7^ 1. 

At next-to-leading order (NLO) there are tree-level contributions from local counterterms 
involving one insertion of the quark mass matrix. In QCD we write the operator set as 



Wn),3 



( bt ] Tr 



B 



+ b<? } Tr I B |{A^ + , B},m qi+ \\ + bf> Tr I 5{{A| + , B}, m g5+ } 



+ 4 n) Tr 

(n) 



B{ 



+ bf ] Tr 
+ b^ Tr 



BB 



Tr 



+ bt ] Tr 



B\\ + B 



Tr[m 95+ ] ) V-^n 



BBX 3 ^ 
- traces 



Tr [ m qiA 



(34) 



where Af + = | (£A 3 £ f + £ f A 3 £) and m gf+ = | (^m^ + ^m^) . The 6j n) must be deter- 
mined, however, isovector charge conservation ensures that bp = 0. The explicit renormal- 
ization scale dependence of the matrix element at each order in the chiral expansion requires 
that, in general, the b^ 1 have logarithmic dependence upon the renormalization scale, jj,. In 
QQCD the analogous operator set takes the form |24| 

Qf)(n),3 \ An) -gkji r y3 , n 

+ 6 W (_)(^)(^-HM) g** { X 3 + j ^ + } n g.. n 

+ 6 W B feii (a 3 + )' (^JBwk 



^3 \» 



nlk 



>nlh 



+ 4 n) 0** (A' + )|By fc str(.M + ) 
+ 4 n) (A^ + )™%„ str (M+) + 4 n) B** B m str (\ 3 + M 
] v^w-V _ traces , 



(35) 



where the 6-™^ of QQCD will, in general, differ from those of QCD. Isovector charge conser- 
vation forces bi = 0. 



At tree-level the matrix elements in an octet baryon "z" is 



(O)tree 



+ M 



(l)tree 



( Q 0^...jT = V m ...v^ ( "M?*™ + Q M 



1 (l)tree 



traces 
— trace 



(36) 



where the Clebsch-Gordan coefficients, Aif ee and ®J\4\ Tee , are given in Tables |, |H| and |TJ. 
Parametrically, but not numerically, the leading order matrix elements in QCD and QQCD 
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Baryon 


^(O)tree ^ Qj^ (O)tree 


P 




£+ 





TABLE I. The leading order tree-level matrix elements of the isovector twist-2 operators 
between octet baryon states. The constants and (5^ are defined in eq. (|3l|). 



B M {1)trcc (QCD) 

~V m (-6i n) + bf ] - bf ] + + + m s (6™ + bf ] + bf ] + bf + bf ] ) 

£+ m (4 n) - 2bf ] + 2b^) + m s (-6^ + b^ ] ) 

~° m (4 n) + 4 n) - bf ] - bf ] - 26< n) ) + m s + 6^ n) + 6< n) - fcf - 6^ n) ) 

TABLE II. NLO tree-level contributions to the matrix elements of O^x^.-nn in the octet 
baryons. m = m u = in the isospin limit. 



B 



Q M (l)tvee (QQCD) 



Tn) jJn) 1 An) 



V2°3 



7^ 
770 



-26 



I m (2^ n) + 5/V 



(«) , li.(n) 
2 + 2°3 

i m (46S n) + 4 n) 



+ 46 2 

(«) ■ 
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■/??■< 
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2b ( r ] + i^ n) ) 



2b 



(n) 



+24 n) ) 



NLO tree- level contributions to the matrix elements of ^O^^...^ 



TABLE III. 

baryons. m = m % 



in the octet 



i n the isospin limit. 



are the same at tree- level. We consider matrix elements in only the proton, S + and H°, as 
matrix elements in the other members of the octet are either trivially related to these, or 
vanish by isospin symmetry. 

In QCD, contributions that are non-analytic in m q first appear at NLO and have the form 
m q \ogm q |J. However, in QQCD there are additional contributions of the form logm 9 to 
baryons other than the nucleon from hairpin interactions. We will consider the two distinct 
contributions separately, and write the matrix element at one-loop level as 

(n),3 \ //s-\(n),3 \trcc , //r\(n),3 \non— HP 

piH2...Hn/i \^H1IJ,2—Hn'i \ Ml /i2 ••■/<•»/ i 

(n)>3 \ _ /Q/n(n),3 \trcc , /QWn),3 \non-HP , iQ/r\(n),3 \HP C\7\ 

HXH2—Hn'i \ /il/i2— /W* ' \ HiH2-IJ.nl i ' \ fllH2 .../in '« ' \ / 

The superscripts on the one-loop contributions indicate whether or not a hairpin interaction 
(HP) appears in the diagram. Of course, there are no HP contributions in QCD. 



1. Non-Hairpin Contributions 

The one-loop diagrams in figs, (fj) give contributions that are non-analytic in m q of the 
form m q \ogm q . The diagrams in figs. ([3]) also give terms of the form m q \ogm q in the limit 
that the decuplet-octet mass splitting, A, vanishes, but in general give a more complicated 
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✓ X^X 
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(c) 



B 
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B 



(d) 



✓ x^.X 

B ' N B 
ft ft 

(e) 

FIG. 2. The meson loop diagrams that give leading non-analytic contributions of the form 
m q \ogm q to the matrix element of^O^^ ...^ n between octet baryon states. The crossed circle 
denotes an insertion of an operator from eq. f fj^ J or eq. (^3^), arising directly from the twist- 2 
operator. The smaller solid circle denotes an insertion of the strong two-pion-nucleon interaction 



from the nucleon kinetic energy term in eq. (\22l), while the square denotes an insertion of the 



axial-vector interaction oc F,D,~f. The label "it" denotes an octet meson, and "x" denotes an octet 
fermionic meson. Diagrams (a)-(d) are vertex corrections while diagram (e) denotes wavefunction 
renormalization. 




(a) 




(b) 




(c) 

FIG. 3. One-loop diagrams with %jk intermediate states that contribute to the matrix elements 
of ^O^/ll .,.fj, n between octet baryon states. The thick solid line inside the loop denotes a %jk 
propagator, while the dashed line denotes a meson propagator. The crossed circle denotes an 
insertion of an operator from eq. Mi 



or eq. ( 31_ ), arising directly from the twist- 2 operator, and 
the square denotes an insertion of the strong interaction vertex oc C. Diagrams (a) and (b) are 
vertex corrections while diagram (c) denotes wavefunction renormalization. 



11 



expression. We write the non-hairpin one-loop contribution to the forward matrix element 
in terms of the contributions from the meson mass eigenstates, 



(0(n),3 v 
(Q(T)(n),3 \ 



non-HP 



non-HP 



1-S nl 
8vr 2 / 2 

1 -5 nl 
8tt 2 / 2 



E M? Hj + £ A/? J,- 

E Q ^? } + E Q K U) * 



— traces (38) 



— traces 



As we are assuming isospin symmetry in our calculation, the sum over mesons for QCD 
corresponds to j = n, K, r], while for QQCD the sum corresponds to j = n, K, r] s . The 
coefficients resulting from mesons with the mass of the pion, Ai^ and ^Ai^ [ Af^ and 
QAfW ], from the diagrams in fig. (^) [ fig. (^) ] are tabulated in Table [TV] [ Table [VII] ], 
while the coefficients resulting from mesons with the mass of the kaon, Av~ K > and ®A4( K ^ [ 
A/W and 1, are tabulated in Table |V| [ Table |V11J| . The pionic contribution to the 



B 



(QCD) 



Q M^ (QQCD) 



P 



- [1 + 3(D + F) 2 ] M^ ee 2D(D - 3F)M% ee - (a^ + pW) 2D(D - F) 
[1 + D 2 + 3F 2 ] M^f - (qW - 2/3< n >) DF D 2 M%f - \ ( a W + 2/?( n )) D 2 

- [1 + 3(D - F) 2 ] Mgg* 



TABLE IV. Contributions to the matrix elements of O 
baryons that are proportional to H^. 



(n),3 
jLtl/i2.../x r 



and ^O^/ia-Mn i n the octet 



QCD matrix element in the nucleon agrees with the calculations of Refs. ||[7|], performed in 
SU{2) L ® SU{2) R chiral perturbation theory, when D + F — > g&. The coefficients resulting 



B 




MW (QCD) 


(QQCD) 


P 


l 

2 


1 + 5D 2 + 3F 2 ] M^ ee + (qW + D(D - F) 





£+ 




\ [1 + 3£> 2 + 3F 2 ] M%T ~ | (« (n) - 2/? (n) ) DF D(D 


- 6F)M%¥ + 4 (c/ n ) + p^) DF 



-\ [1 + 5D 2 + 3F 2 ] M^ e + la^D(D + F) D{D - 6F)M^ C + 2 (a^ + DF 



+\a^D 2 



TABLE V. Contributions to the matrix elements of O^^...^ and ^O^a-.-Mn in the octet 
baryons that are proportional to Hk- 

from mesons with the mass of the i], AiS v ' [ Af^ ], in QCD and the mass of the r] s , ^Ai^ 
[ QAf iVs) ], in QQCD from the diagrams in fig. (|) [ fig. (|) ] are tabulated in Table |V] [ 
Table [TX| ] . The integrals Hj and Jj are given by 



Jj = (m 2 j - 2A 2 ) log 



Hj =m 2 \og[^\ , 



m]\ , nA A7^— ^,_„/ A - V A2 - m K^ 



2A JA 2 - m 2 log 



A + JA 2 - m 2 + ie 



(39) 
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(QQCD) 


P 

£+ 
CO 
















_i ( a M _ 2 /3(»)) D 2 


TABLE VI. Contributions to the matrix elements of oj^^— Mn 
and the matrix elements of ^O^^...^ n proportional to H r)s in QQCD 


proportional to -ff^ in QCD 
in the octet baryons. 


B 




AAW (QCD) 






Q AAW (QQCD) 


P 

£+ 
H° 


-C 2 
-C 2 

-c 2 


27Wf e + f ( 7 (n) -4^) 
" lM-f + I - eg) 
- A ( 7 W - 


-C 2 
-C 2 


^tree + 2 ^(n) _ j 

$M%? + J ( 7 < n > - 



TABLE VII. 

baryons that are 


Contributions to the matrix elements of O^^-.-^n 
proportional to J n . 


and Q 0^™|l2---Mn in tne octet 


B 




ATW (QCD) 






Qj\f( K ) (QQCD) 


P 

£+ 


-C 2 
-C 2 

-c 2 


^H ree + 1 (7 (n) - ^f 1 ) 

§A<§? + l i " 
fA^o ee + f (7 (n) "4^) 


-c 2 " 
-c 2 




5 wtree , 13 / (n) _ £«) " 
6 ;Vl S+ + 18 [J 3 ) 

■ |^|ee + 1 ( 7 (n) _ ^ ) " 


TABLE VIII. 

baryons that are 


Contributions to the matrix elements of O^^...^ 
proportional to Jk- 


and ^O^/la.-Mn m the octet 


B 




AfW (QCD) 






9 A/" (r?s) (QQCD) 


P 

£+ 
E° 


-C 2 
-C 2 




" + 1 ( 7 w - ) 


-c 2 " 







TABLE IX. Contributions to the matrix elements of 0^xi«a— /*n proportional to in QCD and 
the matrix elements of ^O^^...^ n proportional to J Va in QQCD in the octet baryons. 



where A is the decuplet-octet mass splitting and rrij is the mass of the meson, and we 
have retained only those terms that are non-analytic in rrij or are required for the correct 
decoupling in the A — > oo limit. It is clear that for A = 0, Jj = Hj. 

2. Hairpin Contributions 

There are contributions to some matrix elements from one-loop diagrams involving the 
hairpin interactions in QQCD. Explicit computation of the diagrams in fig. |] yields 
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r|,r|^ N r|,r| 
B 1 n_ B 

(a) 



B ' > B 

(b) 



FIG. 4. One-loop hairpin diagrams with %jk intermediate states that contribute to the matrix 
elements o/^O^™)]^ ,..^ n between octet baryon states. The thick solid line inside the loop denotes a 
%jk propagator, while the dashed line denotes an rj = T] Uj( i,a or fj = fj U) d )S propagator. The crossed 



circle denotes an insertion of an operator from eq. (31), arising directly from the twist- 2 operator, 
and the square denotes an insertion of the strong interaction vertex oc C. The cross on the meson 
propagator denotes an insertion of a hairpin interaction with coefficient Mq or a$. Diagram (a) 
is a vertex correction while diagram (b) denotes wavefunction renormalization. 



(Q (n),3 \HP 



1 



2C 2 



a$ ( m~ s 1 I 



16tt 2 / 2 9 

2 / T AA T AA 



-A AN 

us J 



+ m l i 1 



-AA -yAA 

uu us 



traces 



(40) 



where a$ and M are defined in eq. (|]) and C is the axial coupling between the decuplet 
and octet. The integrals I ab are defined to be |L4| I qq , = I{m r]q ,m 71 ,, A, A,/x), where 



J(m 1 ,m 2 , Ai, A 2 ,/i) = == , (41) 



[A 1 -A 2 ][mf-m|] 



with 

r(m,A,/i) 



m 2 - -A 2 " 


A log 


fm 2 \ 


2 


A 2 - m 2 ' 


3 . 






+ 3 





log 



'A - \/A 2 - m 2 + ie' 



A + VA 2 - m 2 + ze 



• (42) 



The Clebsch-Gordan coefficients, hi, appearing in eq. ( flTf ) are given in Table [X|. It is 
important to note that there is no contribution from one-loop hairpin diagrams with Bijk 
(containing the octet) intermediate states, since the vertex diagrams are exactly canceled by 
wavefunction renormalization. Further, the hairpin diagrams do not contribute to the matrix 



B 



hi 



P 

£+ 

WO 



TABLE X. Contributions to the matrix elements of the isovector twist-2 operators 

lA 

from one-loop hairpin diagrams involving 7^ intermediate states. 



3 



elements in the nucleon in the limit of exact isospin symmetry, as the possible intermediate 
states in %j k are inaccessible. 
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C. Matrix Elements in the Nucleon 



Our motivation for undertaking this work is to construct a framework in which experi- 
mental data can be compared with the predictions of lattice QCD. As data exists only for 
the nucleon, it is natural to ask why we bothered with the other members of the baryon 
octet. As the present values of rn 1 ^ 1 ' are large, and in fact m^ att ' ~ m^d , we consider three 
light quarks as opposed to just two. Thus the quantities a^ n \ fi( n \ and combina- 

tions of the b\ n ^ must all be determined, in addition to the strong interaction couplings F, D 
and C in order to perform the quark mass extrapolation at one-loop order. For a one-loop 
analysis as we have presented here, the coefficients 7^"- ) and can be determined at tree- 
level from calculations of matrix elements in the baryon decuplet, and it is consistent to use 
these values in the one-loop expression for the octet baryons. Calculation of the proton, 
S + and 5° matrix elements will allow for a determination of a^ n \ /3^ n \ and combinations 
of the b\ n \ and hence the m g -dependence of the nucleon matrix element. While the hairpin 
interactions, Mq and do not contribute directly to the matrix elements in the nucleon, 
they do contribute to the matrix elements in the other octet baryons, and so will indirectly 
influence the calculation of the matrix elements in the nucleon. 

To get a feel for the m^-dependence of the matrix elements in the proton, we pick some 
values for the parameters describing the proton matrix element and vary the pion mass. 
For the QCD matrix element shown in fig. [|, we choose 7^"- ) = a^ n > = 0, and cr n ) = +2 
, = +1 and b\ n ^ = for n / 1 (this does not eliminate contributions from the 7^ 
intermediate states, as they contribute in wavefunction renormalization) . Also, we use the 
tree-level QCD values for the axial coupling constants F = 0.5, D = 0.8 [^5| and C = 1.8 p6| . 



and the measured value of the decuplet-octet mass splitting, A = 293 MeV. For the QQCD 
matrix elements shown in fig. |^, we choose three different parameter sets, set (a), set (b) 
and set (c), that reproduce the QCD matrix element at m n = 400 MeV. Set (a) has 
ry( n ) = (j( n ) = o, a( n ) = +2.35, (3^ = +1 and = for the tree-level matrix elements, with 
F = 0.5, D = 0.8, C = 1.8, and A = 293 MeV. Set (b) has 7W = = 0, = +2.11, 
fi( n ) = +1 and of" = for the tree-level matrix elements, with F = 0.5, D = 0.8, C = 1.8, 
and A = 500 MeV Q. In contrast, set (c) has = = 0, qW = +0.89, (3^ = +1 
and b\ n ^ = for the tree-level matrix elements, with F = 2.0, D = 2.0, C = 1.0, and 
A = 500 MeV. 

IV. CONCLUSIONS 

It is important to know the parton distribution functions of the nucleon in QQCD. 
In order to obtain these predictions from the lattice computations of the foreseeable future, 
extrapolations in m q from the lattice values down to the physical values will still be required. 
We have taken the first step toward establishing the theoretical m^-dependence by computing 
the leading m^-dependence of the forward matrix element of the isovector twist-2 operators 
in the octet baryons in QQCD using Q%PT. It is certainly true that higher order calculations 



Analysis of quenched simulations suggests that the mass splitting A is ~ 400 - 500 MeV pi. 
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FIG. 5. The variation of the matrix element of the isovector twist- 2 operator in the proton. 
The kaon mass is fixed to its experimental value. The solid curve corresponds to a QCD matrix 
element with = = 0, = +2, (5^ = +1, b\ n) = 0, F = 0.5, D = 0.8, C = 1.8 
and A = 293 MeV. The dot- dashed curve corresponds to the QQCD matrix element with the 
parameters of set (a) described in the text, 7^) = = 0, = +2.35, (3^ = +1, = 0, 
F = 0.5, D = 0.8, C = 1.8, and A = 293 MeV. The dotted curve corresponds to the QQCD 
matrix element with the parameters of set (b) described in the text, 

7 (n) = a (n) = Qj a (n) = +2.11, 

= +1 ; b\ n) = 0, F = 0.5, D = 0.8, C = 1.8, and A = 500 MeV. Finally, the dashed curve 

corresponds to the QQCD matrix element with the parameters of set (c) described in the text, 
7 (n) _ a {n) = 0; a (n) = +Q gg? p(n) = +1> M = Q p _ 2 ^ p _ 2 ^ Q _ ^ and A _ 5QQ Me y 



are required, not only to establish the convergence of the expansion about the chiral limit, 
but also to allow an extrapolation to the relatively large values of the lattice quark masses. 
Even in QCD, the convergence of SU(3) chiral perturbation theory varies from observable 
to observable due to the large value of the kaon mass. Recently, it has been observed that 
baryon masses obtained from QQCD simulations by applying the QCD extrapolation agree 
reasonably well with QCD simulations of the masses [0. It will be a pleasant surprise if 



this empirical feature is found to hold for other observables, and particularly the matrix 
elements of the twist-2 operators. We look forward to a complete analysis of the quenched 
lattice data, including an extraction of the constants that describe these matrix elements 
and eventual predictions at the physical value of m n . 



MJS is supported in part by the U.S. Dept. of Energy under Grant No. DE-FG03- 
97ER4014 and JWC is supported in part by the U.S. Dept. of Energy under Grant No. 
DE-FG02-93ER-40762. 



16 



REFERENCES 



[1] G. Martinelli and C.T. Sachrajda, Phys. Lett. B 196, 184 (1987); Nucl. Phys. B306, 
865 (1988); M. Gockeler et al, Phys. Rev. D 53, 2317 (1996); M. Gockeler et al, Nucl. 
Phys. Proc. Suppl. 53, 81 (1997); C. Best et al, |hep-ph/970650l . 

[2] D. Dolgov et al, Nucl. Phys. Proc. Suppl. 94, 303 (2001); D. Dolgov, Ph.D. thesis, MIT, 
Sep. 2000. 

[3] J. W. Negele, Nucl. Phys. A699, 18 (2002). 

[4] S. R. Sharpe and N. Shoresh, |hep-lat/0108003| ; [hep-lat/0011089| ; Phys. Rev. D62, 
094503 (2000); Nucl. Phys. Proc. Suppl. 83, 968 (2000); S. R. Sharpe, Phys. Rev. D56, 
7052 (1997); C. W. Bernard and M. F. L. Golterman, Phys. Rev. D49, 486 (1994). 

[5] A. W. Thomas, W. Melnitchouk and F. M. Steffens, Phys. Rev. Lett. 85, 2892 (2000); 
W. Detmold, W. Melnitchouk, J. W. Negele, D. B. Renner and A. W. Thomas, 
hep-lat/0103006| ; W. Detmold, D. B. Leinweber, W. Melnitchouk, A. W. Thomas, 
and S. V. Wright, |nucl-th/0104043| . 

[6] D. Arndt and M.J. Savage, Nucl. Phys. A697, 429 (2002). 

[7] J-.W. Chen and X. Ji, [hep-ph/ 01 05197 . 



[8] J-.W. Chen and X. Ji, [hep-ph/ 01 05296 



[9] J-.W. Chen and X. Ji, Phys. Rev. Lett. 87, 152002 (2001). 
[10] S. R. Sharpe, Nucl. Phys. B17 (Proc. Suppl.), 146 (1990). 
[11] S. R. Sharpe, Phys. Rev. D46, 3146 (1992). 

[12] C. Bernard and M. F. L. Golterman, Phys. Rev. D46, 853 (1992). 

[13] J. N. Labrenz and S. R. Sharpe, Phys. Rev. D54, 4595 (1996). 

[14] M. J. Savage, |nucl-th/0107038| . 

[15] S. R. Sharpe, Phys. Rev. D41, 3233 (1990); 

[16] G. Colangelo and E Pallante, Nucl. Phys. B520, 433 (1998). 

[17] M. F. L. Golterman, Acta Phys. Polon. B25, 1731 (1994). 

[18] E. Jenkins and A. V. Manohar, Phys. Lett. B255, 558 (1991). 

[19] E. Jenkins and A. V. Manohar, Phys. Lett. B259, 353 (1991). 

[20] E. Jenkins, Nucl. Phys. B368, 190 (1992). 

[21] E. Jenkins and A. V. Manohar, Baryon Chiral Perturbation Theory, talks presented at 
the workshop on Effective Field Theories of the Standard Model, Dobogoko, Hungary 
(1991); 

[22] For a recent review see U.-G. Meifiner, Essay for the Festschrift in honor of Boris Ioffe, 
in "Encyclopedia of Analytic QCD", edited by M. Shifman, to be published by World 
Scientific. [hep-ph/0007092| . 

[23] M. F. L. Golterman and E. Pallante, JEEP 0110, 037 (2001); |hep-lat/0108029| ; 
[hep-lat/01 10183] ; 

[24] J.-W. Chen and M. J. Savage, |hep-lat/011105C . 

[25] R. Flores-Mendieta, E. Jenkins and A. V. Manohar, Phys. Rev. D58, 094028 (1998). 

[26] M. N. Butler, M. J. Savage and R. P. Springer, Nucl. Phys. B399, 69 (1993). 

[27] R. D. Young, D. B. Leinweber, A. W. Thomas and S. V. Wright, [hep-latt/01 11041. 



17 



